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ABSTRACT 



A frequency dependent analysis of a shielded edge-guided 
mode isolator is presented. A Fourier transform technique is 
applied to the boundary expressions of a structure built on 
a dielectric substrate, and the resulting equations are 
solved for the wavelength ratio. By using perturbation 
analysis and the results obtained for the dielectric case, 
solutions for the normalized propagation constant and 
attenuation for waves traveling in the -Z and +Z directions, 
in a structure built on a ferri te substrate, are obtained. 
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I . INTRODUCTION 



The analysis of various configurations of microstrip 
transmission lines, is of great importance, primarily, to 
the industrial area, since these kinds of transmission lines 
are easy to manufacture, and they are suitable for use in 
microwave integrated circuits due to their small dimensions. 

Two major groups of microstrip transmission lines can 
be defined, related to the substrate material that is used. 
These are the dielectric and ferrite substrate groups. This 
study presents the analysis of a new type of ferrite built, 
shielded microstrip isolator, using the edge-guided mode of 
propaga ti on . 

It was shown by M. E. Hines [Ref. 1] that an edge-guided 
mode of propagation occurs in a wide microstrip transmission 
line using a ferrite slab magnetized perpendicular to the 
ground plane. The R. F. fields patterns in such a structure 
are shown in Figure 1 [Ref. 1]. 

Propagation occurs in both +_Z directions with equal 
phase velocities and loss, and the R. F. fields patterns are 
mirror images for both directions of propagation. In the 
dominant mode the energy is shifted from one side to the 
other with reversal of the direction of propagation. Hines 
suggested that this phenomenon could be used in the analysis 
of nonreciprocal devices such as isolators, phase shifters 
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X 



Figure 1. R.F. Fields of Dominant Mode in Micro- 
strip Line. D.C. Magnetic Field in 
-y Di re ct i on . 
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and other components, where non re ci procal behavior can be 
attained by perturbing the structure in an asymmetrical way. 

To achieve isolation, a resistive load can be placed 
only on one side of the structure; thus high transmission 
loss occurs for the propagation direction for which the 
energy is concentrated on the lossy side, and lower loss 
results for the opposite direction where the energy is 
concentrated on the other side. The disadvantage of this 
field displacement isolator was the small ratio between 
backward and forward losses, which was not sufficient for 
isolation purposes. 

An improved edge-guided mode isolator was suggested by 
K. Araki , T. Koyama and Y. Naito [Ref. 2] for which no lossy 
electric material was used, but instead one edge of the 
conducting strip was shorted to ground. 

This isolation structure, shown in Figure 2, was built 
and tested and the experimental result showed large attenua- 
tion in the backward direction and small insertion loss in 
the forward direction. 

The following study is an investigation of the theo- 
retical behavior of a model for this isolator structure. 

The exact configuration of the isolator model is shown 
in Figure 3. Since most of the field is confined in the 
dielectric substrate then the left, right and upper conducting 
walls of the shield have negligible affect on the electric 
and magnetic fields configuration, thus the model which is 
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Figure 2 . Isolator Structure. 

easier to analyze can replace the open boundary -isolator 
configuration [Figure 2]. The shield is a rectangular wave- 
guide which implies that the highest frequency of operation 
should not exceed the cut-off frequency of the TE-jq mode of 
the waveguide. 

The mathematical method which was selected for this 
study was used in earlier works like references 3, 4, and 
5. This mathematical method can be used for calculating 
wavelength ratio and characteri s ti c impedance for non TEM, 

TE or TM transmission lines. Essentially, this method deals 
with the boundary conditions of a structure, after they were 
transformed into the spectral domain. Then by using the 
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method of moment [Ref. 4] and assuming either current density 
distributions or electric field distributions (depending on 
the structure) numerical solutions for the wavelength ratio 
and the characteristic impedance can be obtained. 

In essence, the analysis procedure is a complete 
solution for the dispersion characteristics of a structure 
built on a dielectric substrate, and then by applying 
perturbation theory, the phase and the attenuation constants 
of a structure built on a ferrite substrate are obtained. 

This ferrite built device is the isolator. 
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Figure 3. Three dimensional Isolator and Regions. 
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II . DISPERSION CHARACTERISTICS ON DIELECTRIC SUBSTRATE 



A. FIELD AND BOUNDARY CONDITIONS 

Assume waves are propagating along the structure shown 
in Figure 3, in the -Z direction such that the propagation 
term is e Jojt+YZ where u is the radian frequency, and y=a+j3 
is the propagation constant. We assume from now on that for 
the dielectric substrate y is pure imaginary, hence ct=0. 

Later we shall see that for the ferrite subststrate, y is a 
complex number including both terms a and 3. 

The longitudinal components of the electric and magnetic 
fields can be expressed in terms of the scalar electric 

p k 

potential <t> and the scalar magnetic potential <p as 
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K Ci 2 


(x,y)e yZ 
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(x,y )e YZ 






(2) 


where Kc^ = K.. 


2 * , 2 


i = 1,2 


defining the 


spatial 


regions of the 


structure (Figure 


3). Using the 


two 


f i e 1 d 


expressions (1) 


and (2) 


together 


with Maxwe 1 1 1 s 


curl 




equations leads 


to the 


trans verse 


components of 


the 


electric 


and magnetic fields as 


f ol 1 ows : 
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( 3 ) 
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e yZ 
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= (y— - jcoe. 
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Tx“^ 


e yZ 


Hy i 
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i = 1 


,2. 




Applyi ng 


boundary conditions 


at the 


interface between 


the dielectric 


substrate and 


the 


bottom 


wall of the 


waveguide 


(y = 0), and assuming perfect 


conducting walls, tangential 


electric fields should be zero or 


expl i ci t ly wri ti ng 




Ex 2 


(x.o.z) = 0 








(7) 


Ez 2 


(x ,0 ,z ) = 0. 








(8) 


At the interface between regions 1 


and 2 , (y = D) 


the 


tangential electric field must be 


conti nuous : 




Ex i 


(x,D,z) = Ex 2 


(x,D, 


z) 




(9) 


Ez i 


(x,D,z) = Ez 2 


(x ,D, 


z). 




(10) 



The electric fields at the interface between regions 
1 and 2 exist only in the dielectric portion of the interface 
and can be written as: 
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Ex-j ( x , D , z ) 




on strip 
(x)e yZ 



W<X<A 



( 0 on s tri p 
(e 2 (x)e yZ W<X<A . 



(ID 



( 12 ) 



The tangential magnetic fields at the interface are 
discontinuous by the surface current densities, and assuming 
i nf i ni tess i mal 1y thin strip (can be considered as current 
sheet) the following boundary conditions can be written at 
y = D 

rO W<X<A 

Hx, (x , D ,z ) - Hx ? (x,D,z) = < (13) 

VJz(x)e y on strip 



(0 W<X<A n4 x 

Hz, ( x , D , z } - Hz« (x , D ,z ) = <^ K > 

(jx(x)e y on strip. 

At the interface between region 1 and the upper wall of 
the waveguide (y=H) the tangential electric fields must be 
zero or explicitly writing 



Ex^x.H.z) = 0 


(15) 


E z i (x,H ,z) = 0. 


(16) 



Substituting the field expressions of equations (1) 
through (6) into the boundary condition expressions of 
equations (7) through (16) yields the following equations: 
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3*,(x,0) 



3 X 



jwy2 



;(x,o) 
Ty" 



= 0 



(17) 



Kc 2 2 $ 2 (x,0) = 0 



(18) 



3 4 ,® ( X ,0 ) _ 3*,(X,D) 3* 2 (X,D) 



34,(X,D) 



3X 



Juii 



1 3y 



JulJ 2 W 



Kc-j 2 <|>®(X,D) = Kc 2 2 4 2 (X,D) 



(19) 

( 20 ) 



3 4 f ( X , D ) 

r 5x i«i*l 



3 4 > ^ (X , D ) (0 on s tri p 



3y 



ex(x) W<X<A 



( 21 ) 



Kc 



* ♦*(X,D) = 



0 on strip 
e z (x) W<X<A 



3 4>V ( X ,D ) 3 4 >? (x ,D ) 

1 + j(oel 1 



3 X 



ay 



( 22 ) 



(23) 



- (y 






3 X 



+ j ue2 



9 <j> o ( X , D ) 



sy 



-)= 



ro w<x<a 

(Jz(x) on strip 



K C] 2 ^(X,D) - Kc 2 <j >2 (X , D ) = 



0 W<X<A 
Jx(x) on strip 



3*®(X,H) 



ax 



- jup 1 



3^(X,H) 

3y 



= 0 



Kc ] 2 4 ®(X,H) = 0. 



(24) 

(25) 

(26) 
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B. SPECTRAL DOMAIN TRANSFORM 

6 h 

The scalar potential functions and must satisfy 
Helmholtz's equations in the two spatial regions, thus 

{V xy + Kc i ) = 0 (27) 

(7 xy + Kc i ) *i( x >*) = 0 (28) 

2 2 2 2 2 2 
where Kc^ = y + = K. - £ and vxy is a two-dimensional 

Laplacian operator. 

As was suggested by Itoh and Mittra [Ref. 3] a Fourier 
transform was used in the a domain. The transform is defined 
by: 

00 

F x [* 1 (x.y)] = # j ( a , y ) = / *.(x,y) e jaX dx (29) 

i = 1,2. 

Since the structure is closed boundary, a finite Fourier 
transform must be used instead of the infinite Fourier 
transform. The finite Fourier transform is given by 

a 

F x [* t ( X .y)] « 4^ (any ) = / 4 i e janx dx (30) 

0 

, 2 • 7T • n 

where an = — ^ — 

and 

a <j> . ( x , y ) 

F x [ — L n — ] = - F x C* 1 ( x .y)]. (31) 
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By using equations (30) and (31), a general transform 
of equations (27) and (28) is given by 



3 <t>? ( x ,y ) ' 3*f(x,y) o 

F x [— L - 2 ] + F x £-772 ] + Kc ? F x ^i (*.y)] = 0 (32) 

3 X ** ** 



-1 

3y 



Or explicitly writing: 



3 $-(ar\y) 2 

2 — = Y ^ $^(an,y) 

ay 



(33) 



where y? = o.^ 2 + B 2 - K i 2 i = 1,2. 



Equation (33) should be analyzed for both regions 1 and 
2, and the solutions will be the transforms of the electric 
and magnetic scalar potentials for both regions. 

For region 1 , 

y } Z = o n 2 + B 2 - K 1 2 (34) 

Where K, = - & and B = fl . 

Substituting and B into equation (34) obtains, 

V = “n 2 + ^ - !]• < 35 > 
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A 1 is the effective structure wavelength and is related to 
the free space wavelength X by 



AO 



> A' 



AO 



(36) 



where er g ^ ^ is the effective relative dielectric constant 
of the device. 

Since e r-j = 1 in region 1, equation (36) can be 
rewritten as 



x o> x 1 = — - . (37) 

By using equation (37) together with equation (35) one can 
find that 

Yi 2 > a 2 (38) 



Therefore is always a real quantity. 
For region 2: 



2 

y 2 = 




2 



(39) 



2 

where = o / y 2 £ 2 = To 2 £r 2 * 

By substituting K ^ and B into equation (39) one can obtain 

^ ■ % 2 - t»" 2 T 2 ' • («) 
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From equation (40) it is clear that y 2 can be either 
a real or an imaginary quantity, depending on the value of 
an . 

Y 2 will be imaginary for 

- ^ * (^) 2 ] < - < ^ [ Vr2 , r2 - (^M (41) 



and will be re al for 



- " “ an < - !<P Ur 2 e r 2 - (?) 



2 tt J / Xo x “ < an < 

- 



xo 



^ 2^2 - (?> 



(42) 



The last two equations should be carefully observed while 
solving equation (33). 

For region 1 equation (33) can be written as 
2 



34> 1 ( an ,y ) 2 

« = y $i(an,y) 

ay 

and the solution has the form of 

^ 1 ( ocn ,y ) = A(an ) cosh y-| *y + B(an) sinhY-|*y 



(43) 



(44) 



the 



For region 2, two solutions do exist correspon di ng to 
real and imaginary values of y^. 
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For y£ real , 

$ 2 (an,y) = C(an) cosh y 2 y + D(an) sinh y 2 y (45) 

and for imaginary 

ii n 

$ 2 (<*n,y)= : C(an)cosY 2 y + jD(an) sin y 2 y (46) 

n 

where y 2 = jy 2 

After knowing the solutions for $-j(an,y) and $ 2 (an,y) , 

e e h h 

one can write the transforms of , <f» 2 , <j>^ , $ 2 in both regions 
as follows. 



Region 1 D<y<H 

$>®(an,y) = A e (an ) coshy ^ (y-D ) + B e (an) sinh y^(y-D) 



(47) 



ftj(an.y) = A h ( on ) coshy, (y-D) + B h (an) sinh y^y-D) (48) 



Region 2 



$ 2 Un,y) = 



^ 2 (an 1 y) =< 



0<_y <D 

,C®(an) coshy 2 y + D®(an) sinhy 2 y , (y 2 real) (49) 

e ii e 11 

Cy(an) cosy 2 y + j Dy ( an ) siny 2 y, (y 2 imaginary) 

(50) 

c|j ( an ) coshy 2 y + D^(an) sinh y 2 y, (y 2 real) (51) 

h 11 h 11 

C T (an) co Sy 2 y + j D-j- ( an ) s i ny 2 y, (y 2 imaginary) 

(52) 
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Superscript (e) indicates the electric field case and 
(h) indicates the magnetic field case. 

All coefficients that appear in equations (47) through 
(52) can be determined by the boundary condition expressions. 
The above can be done by taking the Fourier transform of 
equations (17) through (26), substituting the field 
expressions for both cases of real and imaginary, and 

finally solving for the coefficients. 

The following equations are obtained: 

Hyperbolic Case (y ^ - real) 

-janyC^(an) - j wu 2 Y 2 D [j ( an ) = 0 (53) 



Kc 2 2 C^(an) = 0 



(54) 



Kc-| 2 A e (an) = E’ z (an) 



(55) 




(56) 



Kc-| 2 A^(an) - K c 2 [C^(an) c o s h y 2 D + D|j(an)sinhY 2 D] = 



=A x (an ) 



(57) 
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( 58 ) 



-janyA^(an) + j ue -| y -j B e ( a n) - 
h h 

- {-jany [C[j(an)coshY 2 D + D[j ( an ) s i n h y 2 D ] + 

+ j we 2 ^ Y 2^H ^ an ^ s ^ n ^ y 2^ + y ( an ) cos h y 2 D] } = ^ z (an) 

Kc ^ A e ( an ) = Kc 2 [C®(an)coshy 2 D + D®(an)sinhy 2 D] (59) 

- janyA e (an) - juy ^ y ^ B^ 1 ( a ) = (60) 

- jctny[C^(an)cOShy 2 D + D^(an)sinhy 2 D] - 

h h 

- jcou 2 y 2 [C^J ( an ) S i nhy 2 D + ( an ) COS h y 2 D] 



- jany[A e (an)coshy-j(H-D) + B e (an)sinhy-j(H-D)] - (61) 

- jojuiY-j [A^ 1 (an)si nhy-j (H-D) + B^(an)coshy-j(H-D)] - 0 

Kc^[A e (an)cosh Yl (H-D) + B e (a ) s i nh Yl ( H-D ) ] = 0. (62) 



The coefficients for the hyperbolic case are obtained 
as f ol 1 ows : 



C^(an) = 0 

p E.(an) 

A e (an) = -1—* 

Kc^ 



(63) 

(64) 



D[J(an) = 0 



(65) 



B h (an) = - ^ 2 E z [ an) + j ,„„^ g y (an) 



wU-jY -j KC-j 



wy 1 y 1 x 



( 66 ) 
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( 67 ) 



DuUn) = j - 1 ^,(ctn) 

Kc 2 i sinh Y2 D z 



C H {an) = ‘ 



any 



B e (an) = - 



Kc^ coy 2^2 S ^ 
ctghy-j ( H - D ) 



^ Z ( an ) + ^(opgYgSinhYgD^X 



« x (an) 

( 68 ) 



Kc. 



2 E z (an ) 



(69) 



u an • y • ctghy, ( H-D ) ctghy,(H-D) 

A n (an) = !~2 ff z (an) - j - — ? v (an) (70) 



wy -j y -j Kc-j 



wy -j y -j 



When substituting equation (63) through (70) into 
equations (57) and (58) one can form two sets of equations 



F l j l (an ) 3 )ff x (an) + ? * ( an , 6 ) 2? z ( an ) = / x (an) 



F , H ( an , B ) 2? ( an ) + F J H ( an , $ ) E ( an ) = 3 (an) 



(71) 

(72) 



or in matrix form 



F-j H (an , 6 ) 


F 2 ^ ( an , 3 ) 




\ ( an ) 




«7 X ( an ) 


F 3 H (an ,B) 


( an , 3 ) 




\ (an) 




3 (an ) 

l— z 



(73) 



where y = j3. 



The elements of the matrix [F H ] are given as follows: 
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F.j H (an,6) = - j 



Kc 1 2 ctghy 1 (H-D) Kc 2 2 ctghy 2 D 



UU-J y-| 



wy 2 Y 2 



(74) 



F o (an,3) = j' 



an 3ctghy-j (H-D) c* n 3ctghy 2 D 



wyi y 



in 



wy 2 y 2 



F ^ ( an , 3 ) = - F 2 ^ ( an > 3 ) 



(75) 

(76) 



F A H (an,3) = j 



(an3)^ctghy-j (H-D) aie^^ctghy-i (H-D) 



ajy l Y l Kc l 



Kc n 



" (77) 



y 

a)e 2 y 2 ctghy 2 D (an3) ctghy 2 D 



Kc, 



Kc 



2 a3y 2 Y 2 



The same method can be applied to the trigonometric case 
( y 2 imaginary) and the following equations are obtained: 

c®(an) = 0 (78) 

Dj ( an ) = 0 (79) 

Dy ( an ) = - j ! Z 7 (an) (80) 

1 Kc 2 siny 2 "D z 

Cy(an) = s — E (an) - j — E { an) 

Kc 2 ‘ i utJ 2 y2"siny 2 "D ulI 2 Y 2 " S 1# nY 2" ^ 

(81) 



Kc 1 2 A h (an) - Kc 2 2 [ Cj ( an ) cosy 2 "D + j ( an ) s i ny 2 " D] = 

= j x ( an ) ( 82 ) 
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( 83 ) 



-janyA ll (an) + j ue -j y -j B e ( an ) - j any Cy ( an ) C0Sy 2 " D + 



+ " Dj ( an ) COSy 2 "D = J‘ 2 (an) . 

When substituting equations (69), (70) and (78) through 
(81) into equations (82) and (83) one can obtain t-he matrix 
form 



F j ( an , 6 ) F 2 ( an , 3 ) 
F^ ( an , 3 ) F^ (an , 3 ) 





' ff x (an) 






J 


S' z ( an ) 




J 2 (an) 



(84) 



where again y = j3 

and the elements of the matrix [F T ] are given as follows: 



F ] (an, 3) = - j 



F^ ( an , 3 ) = j 



K Cl 2 ctgh Yl (H-D) Kc 2 2 ctgY 2 "D 



wy -j y 1 wy 2 Y 2 

an3ctghy ^ ( H-D ) an3ctgy 2 "D 



wy- j y -j 

F 2 ^"(an,3) = - F 2 ^"(an,3) 



ajy 2 y 2 



F^ (an, 3) = j 



(an3) 2 ctghy 1 (H-D) ue^^ctghy^H-D) 



u y -j y 1 Kc .j ' 



Kc 



(85) 

( 86 ) 

(87) 

( 88 ) 



(an3 ) ctgy 2 "D °° £ 2 y 2 "ctgy 2 "D 
Kc 2 2 (jp 2 y 2 " Kc 2 2 
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Expressions (73) and (84) can be written in a general 
form as 



F^* T (an,6) Fg ,T ( an , 8 ) 
F 3 ' T ( an > 8 ) • F^ ’ T ( an , 8 ) 



■ 




r / c\ 


E x ( on) 




«F x (an) 


E. (an ) 




J- ( an ) 


L z 




L z 'j 



(89) 



from which solutions to £ x (an) and £’ z (ctn) in terms of 
and ; z (an) can be obtained as follows: 

^ x ( an ) f 4 H * ' (an ,0 ) - e 7 z (an)F 2 HjT (an ,0) 
a x (an) F 1 H,T (an,B)’F 4 H ’ T (an,B) - ,T ( an , 6 ) • F 3 H ,T ( an , 8 ) 

- e7 x ( an ) F^ ( ctn , 0 ) + J z ( an ) • F^ ( an , 0 ) 
r ' ( ’ n| ' F 1 H -' T („ n ,' a) F«. T (.„,.) “ F». r (';„V 6 7^5’ T («r,,e) 

(91 ) 



J x (an) 

(90) 



Define the following terms 



DN - F^ ’ T ( an , 6 ) F 4 ’ T (an, 8 ) - F* ,T ( an , 6 ) F^j ,T ( an , 6 ) (92) 



p H,T 

„H,T _ r 4 ( an , 6 ) 

M 1 ( an , 8 ) * DN 


(93) 


r H,T 

,.H ,T _ ' r 2 ( an , 8 ) 

M 2(an,6) ' DN 


(94) 


F H > T 

M H,T _ 3(an,0) _ mH,T/-\ 

M 3(an,8) ' DN M 2 (on ’ 6) 


(95) 


T . P-] H ’ T (“ n , 8 ) 

4 ( an , 8 ) DN 


(96) 
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which lead to the final expressions 



£ x (an) = ’ T ( an , & ) j an) + M^Un^) j z ( an) (97) 

£ z (an) = Mg *"*"(an , 8) <Z x (an) + M^ ,T (an,B) j z ( an) . (98) 

C. DETERMINENTAL EQUATION 

Nothing seems to be gained so far since neither the 
electric field and the current density components nor their 
transforms are known. In order to simplify equation (97) and 
(98), the method of moments [Ref. 4] is applied in the 
spectral domain. Define a scalar product over the domain 
-®<an<» according to reference 6 as 

n = <» 

<a(an) ,b(an)> = ^ a(an)*b (an). (99) 

n = -co 

Thus when applying this concept to equation (95) and (96) 
and choosing a(an) = £' x (an) or j-^an) and b(an) = «/ x (an) or 
J z { an) respectively, one can obtain 

( an , 6 ) <7 x (an), « 7 x (an)> + ,T ( an , 6 ) j z ( an ) , j x ( an ) > = 0 

( 100 ) 

<M 2 ,T (an ,B)c7 x (an) se 7 z (an)> + <M^j ,T ( an , 3 )j z ( an ) ,j z ( an ) > = 0. 

( 101 ) 
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By using Parseval's theorem it can easily be shown that 



the right hand side of the above two equations is equal to 
zero, due to the orthogonality of £’ z (an) and ^(an), and 
^ x (otn) and ^(an). 

In general, one can expand each current density component 
in a set of basis functions such as: 



In the following analysis, one term approximation is used 
since it was found to be less complicated, without much 
degradation in the accuracy of the final results. 




( 102 ) 



J x (x) = Z a 2k 9zk< x ) • 



(103) 



k= 1 



So with 



i 



J z iar\) = F{al*f z (x)> = a-|f z (an) 



(104) 



and 



I 



•/ x ( on) * F{a 2 .g x (x)> = a 2 -g x (an) 



(105) 



equations (100) and (101) become 




a 2 Zj m i ’ T ( an » 6 ) I 9 X ( an ) I 2 + a i 2 Mil (an ,g ) 



n = -« n = - °° 



(106) 



f z (“ n )g x ( an ) = 0 
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( 107 ) 



a Y m H ,T (a n .s)g x (an) f*(an) + a -j ^M^j’^an.S)' 

2 „ 3 n = — 

n = -oo 

•|f z M)| 2 = 0. 

From the geometry of the device (Fig. 3) and the loca- 
tion of the coordinate system it is clear that all choices of 
current density distributions, J z (x) and J x ( x ) , are neither 
even nor odd functions. 

Each current density distribution can be expressed as 
a linear function of even and odd functions as follows: 

J z (x) = a l [f ze (x) + f zo (108) 

J x (x) = a 2 [g xe (x) + g xo (x)] (109) 

When taking the Fourier transform of equations (108) 
and (109) then according to reference 7, the following 
transforms are obtained 



-7 z (an) = a 1 [f ze (an) 


+ jf 20 (an)] = a 1 f z (an) 


(no) 


=7 x ( an ) s a 2^ g xe^“ n ^ 


+ J 9 x0 (cin)] = a 2 g x (an). 


(HI) 



In order to solve equations (106) and (107) for the 
non-trivial solution, the determinant of the coefficient 
matrix should be set to zero for all sets of physical para- 
meters, at each frequency of operation. The above can be 
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achieved by finding the appropriate 3 that satisfies this 
requirement. From that 3 the dispersion characteristics of 
the device can be calculated. 

The determi nental equation has the form of 



Y ,T ( an , 6 ) | g x (an ) | 2 

n = -oo 



M^ ,T (an,3) I f z (an) I 2 - 

n = -oo 



00 

2Z M 2 ,T ( an > s ) • f z (“ n ) 

n = -oo 



Mj ,T ( an , 8 ) 

_n = -oo 



9x (an ) f z (anl 



= 0 



( 112 ) 



By investigating equations (93) through (96) one may 
observe that M^(an,s) and M^(an,3) are even functions, while 
M 2 ( an , 3 ) and M 3 (an,3) are odd functions. 

Using this information together with equations (110) and 
(111) leads to the final expression of the determi nental 
equati on 



00 

Y (an ,3 ) 1 9 X ( an ) I 2 


• 


M^ ,T (an , 3 ) | f z (an) | 2 


_ n =-oo 




-n = -°° - 



M^^an^) • B(an) 



= 0 



Ln = - 



(113) 



where 

B ( an ) = - f ze (an) • 9 x0 (an) + f z0 (an) • 9 x e ( ct n ) . (114) 



31 



In addition to the requirement that the determinant of 
the coefficients matrix should be set to zero, one can see 
from equations (106), (113) and (114) that the ratio of the 



two coefficients has the form of 





(115) 



,T (otn ,3) B ( an ) 



n = -oo 



Since the expression in the brackets is a real quantity, 
it is well understood that there is a 90° phase difference 
(in time) between the two current density components. 

D. CURRENT DENSITY COMPONENTS 

After the determi nental equation was simplified as shown 
in equation (113), the two current density components J z ( x ) 
and J x (x) were approximated, thus that the equation could 
have been solved. Various one term approximations were 
investigated and the set of components that was chosen to be 
substituted in the determi nental equation was the one found 
as the best approximation of both current distribution 
components . 

The approximated current density component in the 
z-direction has the form 




on the strip 



el sewhere 



( 116 ) 
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and in the x-di recti on 



( a^ cos on the strip 

j (X) = ). (117) 

I 0 el sewhere 

The two current density components are shown in Figure 4. 




Figure 4. Longitudinal and Transverse 
Current Density Components. 

It is clear that the ratio a a -j = e 20 = 4.85 • 1 0 8 
indicates that most of the current that flows in the z-direction 
is concentrated at the right edge of the strip [Fig. 3]. On 
the other hand the choice of J x (x) indicates that there is no 
current in the x-direction at the right-edge of the strip. 
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Computer output indicated a frequency dependent of the 



ratio a^/ a, as shown in Figure 5. This curve was plotted for 

D =0.125" w = 0.45" £ = 16 and a, = 10' 10 . 

r I 




FREQ (GHz) 



Figure 5. Ratio of the Maximum Current 
Density Components 

From Figure 5 one may learn that the transverse current 
component on the strip increases relative to the longitudinal 
component, as the frequency increases. 
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The Fourier transforms of the current density components 
are gi ven as : 



tJ z (an) 




a 1 • f 2 (an ) = 



(118) 



“*( e^cos anw-1 ) + e^ * an • s i nanw 

20 ^ ? 

(j”) + ( an ) 



+ 3 



20 2020 
-an(e cosanw-1) + — e sinanw 

^~2 2 
(^) + (an) 



^ x (an) 



_ / , TT X ianw _ 

- J a 2 COS 2w ‘ 



e = a 2 .g x (an) = 



(119) 



= a„ 2w C0Sanw I_ sinanw + an 

2 ^ 2- + J 5 

_%> - <■"> (f*) - (an)' J 

where the even and odd components of each transform can be 
identified by refering to equations (110) and (111). 

After substituting equations (118) and (119) into 
equation (113), the determi nental equation was programmed into 
Fortran language and the root 3 that set the determi nenta 1 
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equation equal to zero, was found by an iteration method on 
a digital computer. From the value of B, A'/A could easily 
be calculated. 

Figure 6 shows: 

1. Computed and measured A * / A vs. frequency for 
D =0.125", W =0.45", e r = 16. 

2. Computed A 1 / X vs. frequency for D = 0.025 ", 
w =0.45", e r = 16. 

Both curves are plotted for A = 0.9", H = 0.4". 

It can be seen that there is a good agreement between 
experimental and computed wavelength ratio. The highest 
deviation is on the order of 10 percent which can be attributed 
largely to the experimental apparatus. 

One can learn that the wavelength's ratio is directly 
proportional to the dielectric substrate width', which implies 
that the fields are more confined in the dielectric substrate 
as its width decreases. 

E. AVERAGE POWER FLOW 

For later discussion of the characteristics of the device 
built on ferrite substrate, a perturbation technique is used. 
One of the entries to the perturbation expression is the 
average power flow in the device built on a dielectric 
substrate . 

A general expression for the time average power flow in 
the direction of propagation is given by 
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SHIELO DIMENSIONS— 




0.28- 



0.27 H r 

1 2 



3 




FREQ CGHzO 



Fiqure 6. ComDuted and Measured Dispersion Characteri s ti cs 
for D = 0.025" U = 0.45" and D = 0.125" 

U = 0.45" with e =16. 

r 
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( 120 ) 



Pave = i Re // (ExH )*a_ da 
c S 2 

Since (ExH )*a z = ExHy - EyHx* 

and da = dx dy 

then equation (120) can be rewritten as 

Pave = - i Re JJ (ExHy* - EyHx*) dx dy (121) 

c S 

The minus sign was added in equation (121) to make the 
result positive (recall propagation in the -z direction was 
assumed ) . 

Both the electric and the magnetic fields are known in 

e h 

terms of the scalar potentials <j> and <f> thus substituting 
equations (3) through (6) into equation (121) and using 
Y = jB for the waves traveling in the -z direction, leads to 
the following expression 



P i ave =~ 2 Re 



d <P i . d <p . 

( j 3 g x ■ J 



ay 



ff 3*4 3*i 3*® ) 

+ i s J ( j 6 — + (js ~ + dxdy ) (122) 

1 = 1 ,2 

Subscript i denotes the spatial regions, air and dielectric 
substrate respectively, for which equation (122) must be 
eval ua ted . 



38 



The total average power flow in the device is the sum 
of the spatial power components that flow both in the air and 
in the dielectric substrate as shown in the following equation 

Pave ToT = P l ave + P 2 ave (123) 

Applying Parseval's theorem to equation (122) obtains 



P.ave = - 2 a Re 



£j[ 

n = - oo J \ 



an aj3e .• | (an ,y ) | 



(124) 



- an aj0y.|$ 1 . (an,3)| - ai 3 u ^ | 



3*5 ( an ,y) 



ay 



- u 3 e . 



a^?(an ,y) ( 
5y 



j an 3 



( an >y ) — ay 



3^5 ( an >y) 



h * 3^5( an »y) 



ay 



+ j a n K . 






34>^(an,y) 



3y 



.h* (o "’ y) 

+ ay . 



dy 



The limits of integration depend on each region and are 
H D 

J for region 1 and j' for region 2. 

D 0 

After obtaining the general expression for average power 
flow, one can apply it for both regions. 

Recall equations (47), (48), (64), (66), (69) and (70) 
and have the first two modified, the following equations can 
be obtained 
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sin h y i (y-H) 



(125) 



»®(an,y) = A (an ] 



* -J ( an ,y ) = A (an) cos h Yl (y-H) 



whe re 



A (an) = 



£,(an) 



Kc 1 sin h y i ( D-H ) 



A (an) = - j 



ctn6S’ z (an ) 

Kc, 2 



ff x (an) 



OIU -| Y 1 sin h y -j ( 



After using the modified expressions for $®(an,y) and 
$J(an,y) in region 1, and integrating with respect to 
one can obtain 



p i ave ■ h Re 



z 



cu8 



sin h2v, (H-D) 



2 Y 



(H-D) 



? -e 2 ? -h ) 

an e , | A (an) | + Y , Ui|A (an)| 



+ uB 



sinh2 Yl (H-D) 
-e 2 



+ (H-D) 



-h 



+ Yl e j | A (an ) | 



+ jan 



an y-j | A (an ) | + 

e -h* 



B A (an )A (an) 



- K^A ( an ) A (an) 



sinh2 Yl (H-D) 



(126) 

(127) 

(128) 
■H") 

y , 

(129) 
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In region 2, there are two expressions depending upon whether 
Y 2 is a real or imaginary quantity. As in region 1, the 
integration with respect to y can be done an al y ti cal ly , and 
the following equations are obtained. 

For y 2 real , 

QO 

P 2 ave H = 5A Re 2 

n = -oo 

+ y 2 ^2 I c |j( an ) 1 2 



"si nh2 Y2 D 



- D 



+ tuB 



2 

an Uo 



I Cu (an ) I 2 + 



a) 6 



2 

an £ r 



| D^(an ) I 1 



(130) 



+ Y 2 2e 2 I d h ^ an ' ' 2 



K 2 2 D^*(an) C{j(an ) 



" 


si nh2Y 9 D 


“ 




w — — + D 

Y 2 


+ jan 



6 2 D^(an)C^j (an) 



sinh 2 y 2 D 

and for y 2 imaginary, (y 2 = j y 2 11 ) 

00 

P 2 aVe T = ?A Re X/ 

n = -oo 

s i n 2 y 2 " D 



cu 6 



2 i n e i _l ii 3 | r h 1 2 

an e 2 | Dj | + y 2 u 2 I C T I 



D - 



2y," 



+ a) 6 



an 2 u 2 I C t I + 



,i 2 j n e | 2 

+ y 2 e 2 I ' 



D + 



si n2 y 2 " 0 



2 Y," 



+ an 



2 e h* 2 h e* 
$*DjCj + C jDj 



sin2Y 2 D 



(131) 
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I 



D^(an), cj^(an), Dy(an) and Cj(an) are given in equation (67), 
(68), (80) and (81) respectively 

The total power in region 2 becomes 

P 2 ave = P 2 ave H + P 2 ave T (132)- 

Detailed development of equations (129) through (131) 
is given in appendix A. 

Curves of a-domain power distributions for regions 1 and 
2 are shown in Figures 7 and 8 respectively. The ratio of the 
power carried in each region to the total power carried in 
the device for various parameters, is plotted as a function 
of frequency in Figures 9 and 10. 

Several facts can be studied from the graphs. 

(a) The amount of power carried in each region is 
frequency dependent. 

(b) As frequency increases relatively more power is 
carried through region 2, the dielectric substrate. 

(c) Relatively more power is carried through region 1 
as the permittivity is lowered. 

(d) In the low frequency region, more power is carried 
through region 1 as the substrate thickness increases. 
At high frequencies the power carried in region 1 
approaches the same value, for all values of 
substrate thickness. 

(e) As strip width increases more power is carried 
through region 2, the dielectric substrate, 
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Fiaure 7. Average Power Distribution in Region 1 Vs. 

n (an) for f = 4GHz, D = 0.125", W = 0.45", 
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Figure 8. Averaqe Power Distribution in Region 2 Vs. 

n (an) for f = 4GHz, D = 0.125", W = 0.45", 
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Ill . PERTURBATION ANALYSIS OF EDGE-GUIDED 
MODE ISOLATOR ON FERRITE SUBSTRATE 



A. PERTURBATION EXPRESSION FOR PROPAGATION CONSTANT 

Pe rturbati onal analysis is used to determine solutions 
for perturbed problems slightly changed from other problems 
to which solutions are known. The study that was done up 
to this point obtained a solution for the propagation 
constant of an electro-magnetic wave traveling in a struc- 
ture as shown in Figure 3. The goal of this study as was 
stated in the introduction, was to investigate the character- 
istics of waves traveling in the same structure as shown in 
Figure 3 but built on a ferrite substrate. This goal could 
be achieved by using the perturbation approach considering 
the dielectric case as the unperturbed problem, and the 
ferrite case as the perturbed problem. 

A pe rturbati ona 1 expression for the propagation 
constant due to small changes of a material type in a 
guidance structure is given in reference 8 as 



Y 1 + Y* 



ju !S [ ( e.[Ax p ]•£'•£*) + (u.[Ax m ] -H ' *H*) ] ds 

AS (133) 

// ( E*xH 1 + E ' xii*) -a 2 ds 



where primes denote perturbed quantities. 
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When external magnetic susceptibility is used, and 
assuming that the dielectric constant remains unchanged, 
hence [ A x ^ ] = 0 then unperturbed magnetic and electric fields 
are used instead of the perturbed ones, and equation (133) 
can be rewritten as 



jwu 0 // [ 1 * H • H ds 

i s 

IS ( E*xH + ExH*)-a z ds 



(134) 



In the above equation the following notations are used 



y ' = ct 1 + j 8 1 perturbed propagation constant 



y = a + j 3 
AS 

s 

e 



unperturbed propagation constant 

cross section of waveguide which is 
perturbed by a change of the 
material type 

waveguide cross section 

external magnetic susceptibility tensor 



H and E 



unperturbed magnetic and electric 
fields 



The external D.C. magnetic field is applied perpendicu- 
lar to the direction of propagation as shown in Figure 11. 
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Nott: th* uppor wall 
of wavtguidt was 
romovtd* 





Figure 11. Shielded Isolator Structure. 
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According to reference 1 the lines of the magnetic 
field lie in Z-X plane hence an interaction between the R.F. 
field and the ferri te results when the D.C. magnetic field 
is applied perpendicular to Z-X plane. 

The direction of the applied magnetic field implies no 
R.F. demagnetization in the y direction. Thus the external 
magnetic susceptibility tensor is given by 



[x®] = 



X 0 

A xx 



'xz 



K zz 



X - + y - - + 

A xxaxax A xzaxaz 

+ v - - + v - 

A zxazax A zzazaz 



(135) 



and the integrand in the numerator of equation (133) can be 
explicitly written as 

Cx®]-H-H = [ Xxx H x + X XZ H Z ]H* + [x zx H x + X 22 H 2 ]H z * (136) 



The denominator of equation is simply 4Pave so one can 
rewrite equation (134) as 



Jwy. 



D A 

/ /(x xxl H xl‘ 
0 o xx x 



(137) 



+ x zzl H zl 



+ x xz H z H : +x zx H x H z ) dxd y 



4 Pave 
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B. MAGNETIC SUSCEPTIBILITY TENSOR ELEMENTS 

According to small signal approximation of the equation 
of motion in reference 8 two frequencies can be defined as 



u z = K - N y u m + N z “m> = % ’ % 



“x = K ‘ N y “m + N x = ' «V 



(138) 

(139) 



where 

io 0 = yHoy 
“m = ^4 ttM s 

y = 1.4*g[MHz/oe] =0.879*g[Mrad/oe] 
g = Lande' factor ( * 2) 

Hoy = applied D.C. magnetic field in y direction [oe] 

4ttM 3 = saturation magnetization [Gauss] 

N x ,Ny,N z = Demagnetization factors (in this case 
N x « N z = 0. N y = 1) 



Define the ellipticity of the normal modes of the 
uniform precession [Ref. 8] as: 



(140) 



and Kittel resonance relation for an ellipsoid [Refs. 8 and 9] 
as 



“r = V'V'z 
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The Lan dau -L i f s h i tz damping factor is defined in 
reference 9 as 



a 




( 141 ) 



where aH is the line width of the ferrite. 

In this study the ellipticity e Q is equal to 1 and 
\ = “o ■ > so by taking 

X = X' - jx“ ( 142 ) 

and knowing that [Ref. 8] 



X xz “ ~ x zx 



(143) 



then the susceptibility elements can be expressed as 



x xx " x zz 



/ 2 2 W 22 
“m“r^“r ' + “m^r^ a 

A 



X XX " x ZZ 



w^o) a[o)^ + (jo^O+ci^)] 



- u^(l+a^)] 



, , _ , m 

x xz - - x zx ’ J 



y" = - y = i- 

x xz A zx J 



2(j cow a 

m r 



and 



, r 2 2 / , , 2 n -,2 . . 2 2 2 

A = L^ r -a) (1+a )J + 4o) r 0 ) a 



(144) 

(145) 

(146) 

(147) 

(148) 
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Since X xx = x zz and x xz = " x zx then one can sim P lif y 
the perturbation expression given in equation (137) as 

f ol 1 ows 

l Mx xx (|Hx| 2 ♦ |H Z | 2 ) +X xz (H 2 H y *-H 2 *H x ) 
y • + y* = ja)yo 4-pavi (149) 



C. COMPUTATION OF PERTURBATION EXPRESSION IN THE 
SPECTRAL DOMAIN 

In order to solve the perturbation expression [Eq. 149], 
the magnetic fields, H x and H 2 were substituted in terms of 
the scalar potentials as follows 



★ j ojy D A i 9<f>5(x,y) 

+ Y = 4Pl7e " X xx [(±JS- 
00 I 



3 X 



+ j me 2 



3*|(x,y ) 

3y 



) (150) 



(+js- 



3<f>" (x ,y ) 
Tx 



a <{■ | C x »y ) 



4 .h. 



- j(oe 2 — ) + Kc 2 <t> 2 ( x ,y ) ] + 



+ 2x xz [Kc 2 2 ^(x,y)(^6 



3 4>2 c X »y ) 

9X 



- Jwe, 



3 4>| (x ,y ) 

w 



) ] f dxdy.. 



The upper and lower signs denote waves traveling in negative 
and positive z-directions respectively. 

Apply Parseval's theorem 



+ y* = 



4PaveA 



ii /{ x xx [( 6 2 % 2 + KcJ )h h 2 (an,y)| 2 + 
n = *o ^ 



(151 ) 



+ jBana)e2( < ^2^ an, ^^ 



3*o(on ,y ) 



3 $ 5 ( an ,y ) 



h/ - - r ^3y‘ — - *!!*<■" 
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(151 ) 



2 ^ 
“ £ 2 2 



3$ 



|(an ,y ) 



ay 



] + 



.2*h, 



2 x xz Kc 2 $ 2 ^ an ,y )Cl 6 a n * 2 *^ an *y) - Jule 2 - 3y _ ^ 1 dy - 



34; 



ay 



The integration with respect to y was computed 
analytically for both regions where y 2 is either a real or 
an imaginary quantity. 

Recall equations (49) through (52) 



u (c[j(on) cosh y 9 y 

*2 (“ n >y ) H 2 

(c!f (an) cosy 2 "y 



Y 2 rea l 
Y 2 imaginary 



3 $ 2 ( an »y ) 

ay 



Y 2 D^(an) coshY 2 y 

Y 2 "D®(an) cosY 2 "y 



Y 2 re al 
Y 2 imaginary. 



Substituting into equation (151) for both hyperbolic and 
trigonometric cases, and integrating with respect to y, one 
can obtain: 

For the hyperbolic case 

n= *“ ' _sinh2y 2 D 



(152) 



Jwu \ a i n ii y o u p p p u <- 

+ - SPlTiA Z X XX [ 2 y ? + Wl* + Kc 2 )|C H (an)| + 



+ w =2 y 2 1 u H 



2 - 2 2 | D® ( ctn ) j + jeanuY 2 e 2 (c[] (an)D^*(an) 

- cjf(«n) D«(an)] - 



- j 2 X X z Kc 2 u e 2 y 2 D H (an ) [ 



h sinh2y 2 D 



2 Y ; 



+ D] 
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And similarly for the tri gonome tri c case 



( 153 ) 



» * 

Y + Y 



juy 0 

8Pa veA 




M 

sin 2 y 2 D 
2 V 



+ D][(6 2 an 2 



+ Kc 2 2 ) * I C j ( an ) | 



2 



+ u 2 e 2 ^ 1 I ( an ^1 * 6 anu£ 2 Y 2 (Cy(an)Dy (an) + 

+ Cj* ( an ) D® ( an ) ) ] - 

2 « e * h si n2 Y? ' D ) 

- 2 Xxz Kc 2 ue 2 y 2 Dy (an)Cy(an)[ — ii + D] > . 

D. COMPUTATION OF NORMALIZED PHASE CONSTANT AND ATTENUATION 
The perturbed and unperturbed propagation constants have 
the form of 

y' = o' + j 8 1 (154) 

Y = je (ct = 0 for dielectric case) 

Thus the left hand side of the perturbation expression is 
equal to 

Y ' + Y* = o' + j 3 ' + (-je) = a ' + j ( 6 ' - 6) (155) 

Using equation (155) one can obtain 

. 00 
J w y t 1 > 

«' + H*' - - Wl^TA 2> E + J -I J < 156 > 

n = -« 



55 




(a) 




FREQ (GHz) 



Figure 12. a) Attenuation and 

b) Normalized Prooagation Constant for Waves 
Travelinq in +2 and -Z Directions for: 

AH = 75 0 £„ H 0 c = 1916.3 oe, 

4nMS = 1 200 Ga, g = 1 .99 e r - 15.2 
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a (db/cm) 




Figure 13. a) Attenuation and 

b) Normalized Propagation Constant for Waves 
Traveling in +Z and -Z Directions for: 

AH - 75oe, H n r = 2271. 42oe, 

4nMS = 1 200 rfcr g = 1.99, = 15.2 
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Figure 14. a) Attenuation and 

b) Normalized Propagation Constant for Waves 
Traveling in +Z and -Z Directions for: 

AH = 75 oe , Hp £ = 2628.6 oe , 

4-rrMS = 1 200 Ga, g = 1.99, e r = 15.2 
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or 



oo 



M|1 0 

8Pave • A 




(157) 



n = -oo 



CO 




(158) 



Detailed development of the last two equations is given in 
appendix B. 

Figures 12, 13 and 14 were plotted for three different 
values of external magnetic fields Hoy, hence for three 
different resonant frequencies. 

Several facts can be studied from these plots. 

a. The ratio of backward to forward attennuation is 
very large, as required for an isolator. 

b. Both backward and forward attenuations increase 
as the resonant frequency increases. 

c. The ratio of backward to forward attenuation 
increases as the resonant frequency increases. 

d. The bandwidth of the isolator is determined by 
the ferrite linewidth. 

e. The normalized phase constants in both directions 
are frequency dependent and exceed high values as 
approaching the lower side of the isolator's 
bandwi dth . 
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The region of negative values of normalized phase 
constants was not investigated, thus at present no explana- 
tion is obtained. .Hopefully laboratory measurements will 
indicate whether it was an error due to the inherent 
approximations of the perturbati onal technique, or how large 
are the deviations from the correct results. 
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IV. COMPUTER PROGRAMMING 



The computer program that was developed for this study 
was written in FORTRAN IV language and has three major steps: 

1. Computing wavelength ratio A '/A, for waves traveling 
in a device built on dielectric substrate (Fig. 3). 

2. Computing power flow in a device built on 
dielectric substrate. 

3. Computing normalized phase constant gF/g, and 
attenuation in db/cm for waves traveling in a 
device built on ferrite substrate (Fig. 3). 

These computations were done for both forward 
and backward directions of propagation. 

The computations of the wavelength ratio in step 1 were 
made by finding the root S that solved the de termi nental 
equation (113). First, equation (113) was solved for an 
arbitrary value of g, and then, by the use of Newton-Raphson 
iteration method, g was computed until the change between two 
B's in two successive iterations was in the 6th or 7th digit 
after the decimal point. From the value of the last 6, 

A'/A could easily be computed. 

In step 2 the power flow was computed for both hyperbolic 
and trigonometric cases, depending whether was a real or 
imaginary quantity. The computed g, and the coefficient ratio 
a^/a-j obtained in part 1, were used as entries in part 2. 
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In step 3 both power flow and computed 3 were used as 
entries. This step was executed twice each time for waves 
traveling in both directions. The output of this step could 
be obtained for only one value of the following parameter at 
each run: 



All three steps could be calculated for different 
values of frequencies strip's width and substrate's thickness 
in one run. 

The required input data to the program is given in 
Appendix C. 

The limits for all summations in the program were chosen 
by a trial and error method, since there were no definite 
limits that could be pointed out. 

The upper and lower frequencies of operation are bounded 
as follows. 

Lower frequency is bounded by the strip's width such as 



Ferrite linewidth 
Applied D.C. magnetic field 
Saturation magnetization and 
Lande'-g of the ferrite. 



X Q max 



W 



(159) 




or 




C 



(160) 
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Approximated upper bound was applied according to 
reference 10. This reference discusses the case of a 
rectangular guide with dielectric slab perpendicular to the 
electric field as shown in Figure 15. 





H 

Jl 




D , 


A ► 


J \ 



Figure 15. Waveguide Filled with Dielectric 
Slab Perpendicular to the 
Electri c Fie 1 d . 



The guided wavelength is given as 



X 



g 



y - ^i) 

n er o 




( 161 ) 



The cut-off wavelength of this guide is the one that set 
the denominator equal to zero, so from equation (161) one can 
obtain 



A o mi n 





(162) 



an d 



max 



C 

A o min 



(163) 
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V. CONCLUSIONS 



In this study' a theoretical analysis of a model for a 
shielded edge-guided mode isolator was presented. 

The analysis was based on a complete solution for an 
unperturbed problem (the dielectric case) and then by the 
use of a perturbation technique, a solution was derived for 
a perturbed problem, (the ferrite case). 

Final results indicated that above and below the 
resonant frequency the forward attenuation is very low - 
negligible, while the reverse attenuation is high. These 
are good regions in which to operate as an isolator. One 
result which is still unexplained is that the normalized 
phase constant is negative in the frequency region above 
the resonant frequency. 

A computer program was developed for all steps of the 
theoretical analysis. 
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APPENDIX A 



AVERAGE POWER FLOW IN REGIONS 1 and 2 



Equation (129) indicates the total power flow in region 
1. In order to solve this equation both A e and A h should be 
known explicitly. 

Recall 



A e ( n) 



E z {a n) 

Kc-, 2 sinh Yl (D-H) 



(A-l) 



anBff (an) . 

- *„(«»)] rq—7- r r r- n- r (A-2) 

l^ c ^2 x' J «ni^ Yi si n h y ^ ( D-H ) 

M-| H ’ T ( an , B )«T X (an ) + M 2 H ,T ( an , 6 ) J z ( an ) ( A- 3 ) 

F z (an) = M 3 H,T (an,B)/ x (an) + M 4 H ’ T ( an , B U z ( an ) (A-4) 

Since M^’^(an,B), M 2 ^’^(an,B), ’^ ( an , B ) and 

H T 

* ( an , 6 ) are imaginary quantities, one can define 



A h (an) = 
® x (an) = 



"A 


,T (an . 


to 

II 


« " 
n 2 


,T ( an , 


II 

cn 


M H 


,T (an : 


,0) = 



jm 1 H,T (an ,B) 
jm 2 H ’^ ( an , B ) 
j m 3 H ’ ^ ( an , 6 ) 



( A - 5 ) 



M 4 H,T (an,B) = jm 4 H,T (an ,6) 
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Furthermore, referring to equation (115) one can assume 
that a 2 is a real quantity and a -j is imaginary quantity. 
So one can write „ 



a-j = ja 1 ( A-6 ) 

where a^ is real quantity. 

Substituting equation (A-5), (A-6), (110) and (111) 
into equations (A-3) and (A-4) one can obtain 



£ x (an) = E (otn) + jtf^an) 

^ (an) = E (an) + j E (an) 
z Z 1 1 2 



where 



(on) 

X 1 



E (an) 
x 2 



E. (an) 
Z 1 



a i [ - — m^ ,T (an ,3)g XQ (an ) - 
a l 

m 2 H ,T (an ,6 ) f ze (an ) ] 

• m 1 H,T (an,3) g x0 (an) - 
a l 

m 2 H,T (an ,3) f ZQ (an)] 

a ^ [ - m 3 H ,T (an ,3) g xQ (an) “ 
a l 

m 4 H,T (om,3) f ze (on)] 



(A-7) 
( A - 8 ) 



( A - 9 ) 



( A- 1 0 ) 



( A- 11 ) 
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( A — 1 2 ) 



E. (an) 
z 2 



a-|[~ m 3 H ’ T ( an >s)g xe (an) - 
a l 

m 4 H .’ T (an ,e)f 20 (an)] 



and 

99 

^m 2 H>T (an,6)B(an) 
a 2 _ n= -°° 

a, A u T 2 

,T (an ,6) | g x (an) I 

n=-» 



(A-13) 



After investigating equations ( A - 9 ) through ( A - 1 2 ) it 

could be seen that E (an) and E_ (an) are even functions, 
X 1 Z 2 

while E (an) and E (an) were odd functions. Substitution 
• x 2 Z 1 

of these functions into equations ( A - 1 ) and (A-2) led to 
explicit expresions of A (an) and A (an), hence the power 
flow in region 1 could be calculated. 

It is clear that the power flow in region 2 for both 

the hyperbolic and trigonometric cases could be calculated 

6 h 

in the same way since the coefficients D^(an), Cjj(an), 

o h 

Dy(an) and Cj(an) are functions of E x ( an) and £' z (an) which 
their explicit expressions are given in equations (A-9) 
through ( A - 1 2 ) . 
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APPENDIX B 



NORMALIZED PHASE CONSTANT AND ATTENUATION 



As was stated'in equation (157) and (158) 



uy 0 

8PaveA 






( B - 1 ) 



00 



S ' /S = 1 + 8PaveA R E 



(B-2) 



n = -oo 



So in order to calculate a' and 8'/8, both I and R £ should 
be known explictly. 

For the hyperbolic case, equati on (151) had to be divided 
into real and imaginary components. 

Recall equations (67) and (68) and substrate equations 
( A- 7 ) and (A-8) for £ (an) and # z (an) one can obtain 




' { r _ n 

ajp 2 Y 2 s 1 n l r| Y 2 ^ Kc ^ z 2 



E (an) + £ (an)] + (B-3) 

z 2 x 2 



8 a 

+ j[— V*. (an) - E (an)]} 

Kc 1 1 




lE. (an) + j£ (an)] 
z 1 z 2 



(B-4) 
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Substituting equations (B-3) and (B-4) into equation 
(152) and separating to real and imaginary parts, the 
following equations are obtained 

2 2 2 

r«n £ + Kc/ Q n 2 

RI 1 = ~ 2 2 2 { t“ E 1 («") + E y («")] + ( B - 5 ) 

^ u 2 y 2 K ^ 2 2 



+ 1^2 E 7 {an ) - (on)] > + 

Kc^ Z 1 X 1 



2 2 2 

U) e 2 Yo p p 

4 [# (an) + E (an) + 

K C 9 Z 1 Z 2 



i 2Ban r Ban r t-t 2/ \ , 2/ > n 

+ j { j IE (ctn) + E (an)] - 

KCg K 2 Z 1 z 2 



- E (an )E (an) - E- (an) (an)} 



Z 1 X 1 



R I 



[£ y (an)ff (an) - * (an)r (an)] (B-6) 



2 y 2 L x t ' — z 2 



x 2 



ctghYoD D „ x* x' 

RE = [ + 1 H-* 2 - RI, + -p RI 2 ] (B-7) 

y 2 sinh y 2 D ^ z ^ 



Im = [ 



ctghy 2 D 



][- Hr- RI 



y 2 sinh 2 y 9 D 2 



^ - -f 1 RI 2 ] ( B - 8 ) 
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Similarity for the trigonometric case the coefficients 
h s 

Cy(an) and Dy(an) can be expressed in terms of # x (an) and 
£ z (an) as follows 



c}(an) = 1 {[^V («n) - E (an)] + (B-9) 

a)u 2 Y^sin y p Kc 2 z 2 x 2 



+ j[- San E (on) + s (an)]} 
KCg 2 1 X 1 



Dj (on ) 



j o 11 5* 2 ( an ) 

Kc 2 s i n y 2 D 



( B-l 0 ) 



and both the real and imaginary parts of equation (152) have 
the form of 



R I 



1 



Q 2 2 , vr 2 
8 an + KC 2 

2 “2 “2 
w y 2 




E. (an) 
z 2 



2 

E (an).] + 

x 2 



(B-l 1 ) 



+ [- (an) + E (an)] 2 } + 

Kc 2 ^ 1 1 



2 2 "2 

w £ o Y o 9 9 

+ ^ 2(an) + E 7 + 

KC 2 4 Z 1 z 2 



- 2 _ 3 _i n {Ban r 2^^ + 5- 2 ( a n)] - 

Kc 2 Kc 2 Z 1 z 2 

- £ (an)£ v (an) - E (an)# (an)} 

1 1 X 1 z 2 x 2 
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( B - 1 2 ) 



R I 



2 u 2 x 1 



[ff (an)g (an) - 2 v (an )2 (an)] 



x 2 Z 1 



ctg Y? D x‘ 



-X v . 



r e = [ — 7 -*^ + -- n ] [-r 1 R in + ~-p ri 2 ] (b-13) 

sin i Y 2 D y 2 z \ c . c . 



I B = [' 



si n 



n CtgYp D Xy y Xy 7 



■yJO v 2 



( B-l 4 ) 



After substituting properly equations (B-7), (B-8), (B-l 3 ) 
and ( B - 1 4 ) into equations (B-l) and (B-2) the later can be 
solved for the attenuation and normalized phase constant. 
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APPENDIX C 



COMPUTER PROGRAM 



THEORY OF FIELD DISPLACEMENT DEVICES 

5HEILDED EDGE GUIDED MCCE ISOLATOR 
(BOX'S DIMENSION 0.4' '*0.9'') 



THIS PROGRAM HAS THREE MAJOR STEPS CORRESPONDING TO 
THE STEPS OF THE THEORETICAL ANALYSIS OF THE DEVICE. 
THE THREE STEPS APE: 

1. COMPUTING WAVELENGTH PATIO OF THE DEVICE BUILT ON 
DIELECTRIC SUBSTRATE. 

2. COMPUTING POWER FLOW IN THE DEVICE. 

3. COMPUTING NORMALIZED PROPAGATION CONSTANT AND 
ATTENUATION IN CB PER CM OF FERRITE 3UIL T OEVICE 
BY USING PERTJR BAT I ON THEORY ANC CATA OBTAINED IN 
STEPS 1 £ 2. 



PROGRAM ACCEPTS FOLLOWING DATA: 

1. FIRST CAAD-Ll , M 1, MN 1 FORMAT ( 3 ( I 2? X2 ) ) 

LI -NUMBER OF FREQ. CATA CARDS. 

Ml -NUMBER OF STRIP WIDTHS DATA CARDS. 

MN1 -NUMBER OF SUBSTRATE THICKNESS DATA CARDS. 

2. SECOND CARD- ( FR EQ ( L ) , L= 1 , L 1 ) FG PM/T <D9. 2 ) 

FREQ(L) ARE THE FREQUENCIES IN HERTZ AT WHICH THE 
COMPUTATIONS ARE EXECUTED. 

NUMBER OF ''SECOND'* C ARDS -L 1. 

3. THIRD CARO-(W(M) ,M=1,M1) FORMAT (DIO. 4 ) 

W (M ) ARE THE VARIOUS STRIP WIDTHS, IN METERS. 
NUMBER OF "THIRC" CARCS = M1. 



4. FOURTH CARD-(D(MN),MN=1,MN1) FORMAT (D10.4) 

0 ( MN ) ARE THE VARIOUS SUBSTRATE THICKNESS, IN METERS. 
NUMBER OF "FOURTH'' CARDS =MN 1 • 



5.-IF7H CARD-EPSR ^0 R* AT ( F5 . 2 ) . 

EPSR- RELATIVE DIELECTRIC CONSTANT. 



6. SIXTH CAkD-CELTH FG RM AT ( D12 . 3 ) 

DELTH-L I NEW! DTH OF FERRITE IN OERSTEDS. 



7. SEVENTH CARD-HO >=ORM AT ( 012 . 3 ) 

HO-APPLIED D.C MAGNETIC FIELD IN OERSTEDS. 



8. EIGHTH CARQ-AMAGS FORMAT (D12.3) 
AMAGS-SATURATI ON MAGNET I ZATI ON IN 



GAUSS • 



9. NINTH CARD-G FO RM AT ( F4 • 2 ) 
G- LANDE-G OF FERRITE. 



OUTPUT : 

1. WAVELENGTH RATIO FOR CIELECTRTC SUBSTRATE. 

2. NORMALIZED PROPAGATION CONSTANT FQR BACKWARD AND 
FORWARD DIRECTIONS OF PROPAGATION. 

3. ATTENUATION FOR BACKWARD AND FORWARD DIRECTIONS OF 
PROPAGATION. 



THE OUTPUT IS PRINTED FOR ALL FR E C’JENC I E S , STRIP 
WIDTHS AND SUBSTRATE THICKNESS. 

CAUTTON: THE UPPER ANC LOWER FREQUENCIES OF OPERATION 
ARE BOUNDED TC ELIMINATE UNOESIREC MODES OF 
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propagation in the device. 

ALL COMPUTATIONS ARE PERFORMED IN DOUBLE PRECISION. 
LANGUAGE-FORTRAN 

PROGRAM DEVELOPED BY LT.J.G RAM SFARON ISRAELI NAVY. 
THESIS ADVISOR- PROF. CR. J.3.KNQRR. 

NAVAL POSTGRADUATE SCHOOL MONTEREY CALIFORNIA 93940 
FEBRUARY 1976 

IMPLICIT REAL*8 (A-H.C-Z) 

DIMENSION T0T(4),FREQ(1C),W(8) ,C(8) 

ZZ=20 • DO 
S=0. DO 



FORMAT (D9.3) 
FORMAT (010. 4) 



FGRMAT (3(12, 2X ) ) 
FORMAT ( F 5. 2 ) 

10 FORMAT (C12.3) 



7 FORMAT (F4.2) 

READ (5,3) L 1 » M 1 , MN l 
READ (5,1) (FREQ(L) , L = 1 ,L1 ) 

READ (5,2) ( W ( M ) * M = 1 , M 1 5 
READ (5,2) (D(MN) ,MN=l,Mi\l) 

READ (5,5) E P SR 
READ (5,10) DELTH 
READ (5,10) HO 
READ (5,10) AMAGS 
READ (5,7) G 

P!=3. 1415 926 54D0 

AMUl=4.D0* p I*1.3D-7 

AMUR= 1 • CO 

AMU2=AMU1*AMUR 

EPSl=i.D0/(3fc.D0*PI*1.0+9) 

EPS2=EPS 1*EP$R 

$T EP1-WAVEL ENGTH RATI O CALCULATIONS, 



H-BCX'S HEIGHT 



1 • 0 1 6D-2 



A-BCX'S WIDTH 
A=2. 2860-2 



FREQUENCY LCOP 

CO 4000 L=1 , L 1 
FREQl=FREQ(L)/1.0D+9 
3MEGA=2.C0*PI*FR£G(L) 
0MEGA1= OMEGA/l.OQ+9 
WRITE (6,4) , _ 

4 FORMAT (//////////) 

PRINT OUT QF INPUT DATA 



199 
2C0 
2 C 1 
202 
2C3 
207 



WRITE (6,199) 

PHRMAI <"0X*22H ********** ************> 
WRITE (67 200 ) FREQ 1 

FORMAT (20X,3H* FP EQ1 = » F5 . 3 , 7 H GHZ,2H 
WRITE (6,201) OMEGA 1 

FORMAT (30X,9H* OMcGA 1= , F6 » 3 , 5H GRP--5,2H * 
WRITE (6,202) H 

FORMAT (30X,4H* H=,cll.4,5H MTR , 2H *) 

FORMAT ( ( 30X^4H* A=,S11.4,5H MTR , 2H * ) 
WRITE (6,207) F.PSP 

FORMAT (3QX,7H* E p SR”,F4.1»9X,2H * ) 



* ) 
) 
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WRITE (6,206) DELTH 

206 FORMAT (30X,10H* 0ELTH= ,F7.2,5H CE *) 

WRITE (6,210 HO 

210 FORMAT (30X,10H* C.C • FLD = , F7. 2 , 5H OE *) 

WRITE (6,212) AMAGS 

212 FORMAT (30X,10H* S AT . M AG= , F7 . 2 , 5H GA *) 

WRITE (6,213) G 

213 FORMAT (3CX,10H* LANDE-G= ,F4. 2 , 7X ,1H* ) 

W C I T E (6,208) 

208 FORMAT ( 30X , 22H***** ************ ***** // ) 

SET HIGHER FREQUENCY OF OPERATION 

GM=2*D0*A/DSQRT< l.DO-D(MN)/H*< l.DC-l»DC/EPSR)) 
FREQU=3 • D8/GM 

IF (FREQ(L).GE.FREQU) GC TO 400C 
STRIP WIDTH LOOP 
00 5000 M=1 , M 1 

SET LOWER FREQUENCY OF OPERATION 

FREQL= ( 3 .0C+08/DS GHT ( EPS R ) )/(4*W(N) ) 

IF (FREQ(L).LE. F9EQL) GC TO 5000 
WRITE (6,205) W(M) 

205 FORMAT ( 10X , 2HW = , FI 1 . 5/ / ) 

SUBSTRATE THICKNESS LOOP 



DO 6C00 MN = 1, MN1 
WRITE (6,209) 

2 09 FORMAT ( 20X , 14HLAM8D A 1 /LAMBDA , , 1HC , 15 X , 4H8ETA/ ) 



ni=ioi 

11=201 



BETA-PROPAGATION CONSTANT- 



C 

1GC0 

C 



600 

700 



C 



16 

9 



C 



6 



DB6TA*1.0D-1 
BE T A1 =300*00 

TOT ( 1 >=20. DO 
00 2000 K=1 , 2 



BETA = BET A 1 

IF (K.EC.2) GO TO 600 
GO TO 700 
6ETA=3ETA+08ETA 
AK1=0MEGA*D$QRT ( AMU1*EP$1) 
AK2=QMEGA*DSQRT( AMU2*EPS2) 
AKC1S=AK1**2-BETA**2 
AKC2S=AK2**2-BcTA**2 

SUM 1=0 *000 
SU*M2 = 0 • DOO 
$UM3=Q.DG0 
TTPWR 1 = 0 • DO 
TTP WR 2 = 0 • DO 
TQTP WR=0 • DO 
GO TO 9 



N 1 = 5 1 
11=101 
RF 1 1 =0 • 00 
RE22=0 • 00 
ZM11=0* DO 
ZM21=0 • DO 
N=-N 1 

DO 100 1=1,11 
N=M+1 

ALFAN = 2 .DO*PI*N/ A 

GAMA 1 = 0 SORT (ALFAN**2-AKC1S) 
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GAMA25 = ALFANi**2-AKC2S 
GAMA2=DSQRT (DABS( GAMA2S) ) 

A1 = 1.DG/DTANH(GA,MA1*(H-C<MN) ) ) 

A2=0MEGA*AMU1*GAMA1 

A3=GMEGA*A.MU2*GAMA2 

A4=ALF AN*BET A 

A5=0MEGA*E?S1*GAMA1/AKC1S 

A6=0,MEGA*EPS2*GAMA2/AKC2S 

IF (GAMA2S. GE.O.DCC ) GO TO 20 

TRIGONOMETRIC CASE 

A 7= DT AN ( GAM A2*D ( Mi\ ) ) 

FI =-AKC 1S*A1/A2+AKC2S/ ( A 7* A3) 

F2=A4*A1/A2-A4/ (A 7* A 3 ) 

F3=-F2 

F4=A4**2*Ai/ ( A2*AKC15 ) - A5* Al- 44** 2/ ( AKC 2S* A 7* A 3 ) - 
GO TO 30 

HYPERBOLIC CASE 

20 A8=1.D0/DTANH(GAMA2*D(MN ) ) 

F1=-AKC1S*A1/ A2-AKC2S*A3/A3 

F2=A4*A1/A2+A4*A8/A3 

F3=-F2 

F4=A4**2*A1/ ( A2*AKC1S )-A5*A1+A4**2*AS/< AKC2S*A3 )- 
30 CN=F1*F4-F2*F3 
AM 1 = F4/ ON 
AM2=-F2/0N 
A M 3 = - A M 2 
AM4=F1/DN 

SN=DSIN(ALFAN*W(M) ) 

CN-DC OS (ALFAN+W (M ) ) 

R 1 = DE X P ( ZZ ) 

R2=(ZZ/W(M) )**2+ALFAN**2 

R3=R1*CN 

R7=R1*5N 

Q= ( ( P I / ( 2 • D 0 * W ( M ) ) )**2-ALFAN**2) 

GX 1 = P I / ( 2. D0*W ( M ) ) *CN/C 
GX2=( PI / (2.00*W(M) ) * SN+ ALFAN ) / C 
FZ1=( ZZ/WGM) *(R3-1.D0) +ALFAN*.R7) /P2 
FZ2 = ( -ALFAN* ( R3-1 .DO ) +ZZ/W( M ) *R7 ) /R2 
IF (DABS(TOTd) ) .GE.l.D+01) GO tq 3 

STEP 2 -POWER FLOW CALCULATIONS. 

AMP1=1.D-I0 
Zl = -S UM 2/SUM3 
Z3 = Z 1*G X 1 
Z2=Z1*GX2 

EZ1=AMP1*( -AM3*Z2-AM4*FZ1 ) 
EZ2=AMP1*(AM3*Z3-AM4*FZ2) 

EX1=AMP1*( -Z2*AM1-AM2*FZ1 ) 
EX2=AMP1*(Z3*AM1-AM2*FZ2) 

EZS=EZ1**2+EZ2^*2 
EX S=EX1**2+EX2**2 
IF (S.EQ.1.D0) GO TO 13 

REGION 1 FREE SPACE 

V0=GAMA1*(H-C(MN) ) 

Vl = OMEGA*BcTA*EPSl/2.!:0 
V2=1.D0/(0TANH( VO^GAMAl) 

IP (VO. GE. 85. DO) C-0 TO 12 
V3 = (H-D(MNn/((0SINH<V0) )**2) 

GO TO 11 
12 V 3 = 0 . D 0 

11 V5=BETA/(2.D0*0MEGA*AMU1) 

V4*ALFAN*8ETA/AKC IS 

V 6 -( V4*EZ2-EX2) **2+(-V4*SZl+EXl )**2 

V7=V4* E Z S-E Z 1*EX1-EZ 2- EX2 

TR 1 1 =-V 1*ALFAN**2 /(AKC 15**2) *EZS* ( V2-V3 ) 

TR 12=-V5*ALFAN**2/( GAMA 1**2 )*V 6 *( V2+V3) 
'R13=-V5*V6*( V2-V3) 



A6/A7 



A 6 * A 8 
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14=-V1*(GAMA1**2)/ ( AKC1S**2 ) *F Z3* < V2+ V2 ) 

156=ALFAN/< OMEGA* AM U1*AKC1S ) * ( 3 ETA**2 + AK 1**2 ) *V7*V2 
R1=-(TR1 1+7R12+TR13+TR14+T5 156 ) / ( 2* DO* A ) 

( r.A 11 A 1 C 1C n TO o 



TR 

TR_. 

PWR1 . 

IF (GAMA2S.LE.0 



.00) GG TC 300 



REGION 2 HYPERBOLIC CASE 

V10=GAMA2*0(MN) 

Vll=0MEGA*BETA*EPS2/2 .CO 
VI 2= 1.00/ ( D TANH (V10)*GAMA2) 

IF (VIO.GE.UO.DO) GG TC 14 
V13=0(MN)/( <DSINH(Vi0 ) J**2) 

GG TO 15 

14 V13=0.D0 

15 V15=BETA/(2.D0*GMEGA*AMU2) 

V 14=4 LF AN*6 ETA/ AKC2S 

V16=( V14* = Z2-EX2)**2+(-Vl4*EZl+EX 1)**2 
VI 7= ( V14*EZS-EZ1*EX1-EZ2*EX2 ) 
TR21=-V11*4LFAN**2*EZS/(AKC2S**2) *< V12-V13) 

TR22 = -V15*ALFAN**2/GAMA2S*V16*< V12 + V13 ) 

TR23=-V15*V16*( VI 2- VI 3) 

TR24=-V11*GAMA2S/ (AKC2S**2)*( V12+V13)*EZ$ 

TR256=ALFAN/( OMEGA* AM U2*AKC2S ) * ( 3 E7A**2+AK2**2 ) *V 17*V12 
GO TU 801 

REGIGN 2 TRIGONOMETRIC CASE 

3 00 V25=1.D0/(AKC2$*D$IN(GAMA2*Q(MN) ) ) 
V20=CMEGA*BETA*E3S2/2.DC 
V21=GMEGA*BETA*AMU2/2.0Q 
GET 1=V25*EZ2 
CET 2=- V 2 5*E Z 1 
0ETS=0ET1**2+DET2**2 

V26=ALFAN*BETA*V25/(QMEGA*GAMA2*AMU2) 

V27=1.00/(QMEGA*A.MU2*GAMA2*0S I M G AM A2*D ( MN ) ) ) 

CHT1=-V26*EZ2+V27*£X2 

CHT2=V26*EZ1-V27*EX1 

CHTS=CHT1**2+CHT2**2 

V28=0SI N(2.00*GAMA2*D(MN) ) / ( 2 , DC*GAMA2 ) 

TR21 = -ALFAN**2*V20*0ET$*(D<MN ) -V2 3 ) 
TR22=-ALFAN**2*V21*CH7$*(D(MN ) +V2 8 ) 

TR23=- (GAMA2* 5is 2)*V21*CFTS*(D( M N ) - V 2 8 ) 
TR24=-(GAMA2**2)*V20*CETS*<D( Mi\)+V28) 

TR 25t> = ( BETA**2+AK2**2 )*ALFAN/ 2. CQ*( 0ET1*CHT1+DE72*CHT2 ) 
8 G 1 PWR2=-(TR21+TR22+TR23+TR24+TR256 ) / (2. 00* A ) 

TCTAL POWER FLOW IN FREE SPACE REGION 



C 

C 



C 



TTPWP1=TTPWR1+PWR1 

TCTAL POWER FLOW IN CIELECTRIC PEGICN 
TT P W R 2 = TT P W R 2 + r W R 2 
TOTAL POWER FLOW 

TOTPWR=TCTPWR+PWP 1+PWR2 
GO TO 100 



8 FZS=<R1**2-2.D0*R3+1.CC)/R2 

GXS= ( (PI/(2.D0*W(M) ) )**2+ALFAN**2+ALFAN*oI/W(M)*CSiN 
x< alfan*w(m) ) )/g**2 
T1N=AM1*GXS 
T2N=AM4*FZS 
6=-FZ 1*GX2+F Z2*GX1 
T3N=AM2*3 
SUM1=SUM1+T IN 
SUM2 = SUM2+T 2N 
SUM3= SUM 3+T 3N 
100 CONTINUE 



IF (S. 
IF (TO 



EQ.1.C0) GO ^0 19 
TPWR.NE.0.00 ) GC TC 17 
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SUM3S Q= ( SUM 3 ) ** 2 
TOT ( K ) =SUM1*SUM2+SUM3SC 
It <DABS(T0T(1)I.LE.1.C+C1) GO TO 3000 
. 20 00 CONTINUE 

DT0T=( TOT (2) -TOT ( 1) ) /CEETA 
BETAl=DAeS( BETA 1 -T 0 T ( 1 )/DT3T) 

60 I p (3ETA1.LE.AK1 ) GO TC 50 
GO TO 70 

50 BETA1=BETA1+1.D0 
GO TO 60 
70 GO TO 1000 
30C0 WVLNG=3.0D+8/FREQ(L ) 

RATI 0=2. 00* PI /<3ETA1*WVLNG) 

WRITE (6,211) RATIO, D(MN) ,BET A 
211 FORMAT (24X,F6.4, 7X, E11.5,7X, Ell. 

GO TO 16 

17 S=1 • DO 

SIGN=-1.D0 , PROPAGATION IN 6ACKWARC DIRECTION ( + Z) 

SIGN—1.D0 
8C00 N 1= 1C 1 
11=201 
GC TO 9 
21 S=0 . DO 
6000 CONTINUE 
50C0 CONTINUE 
4000 CONTINUE 

GO TO 8001 

STEP3- CALCULATIONS OF NORMALIZED PROPAGATION CONSTANTS 
ANC ATTENUATIONS. 

18 V14=BETA*ALFAN/AKC2S 

V16=< V14*EZ2-6X2) **2 + ( - V14*EZ 1 + cX 1 ) **2 
VI 7= < V14*EZS-EZ1*EX1-EZ2*EX2) 

V30=( 3ETA*ALFAN ) **2+ AKC 2S**2 
V31=(0MEGA*AMU2*GAMA2) **2 
V22=(CMEGA*EP$2*G AM A 2 ) **2 
V33=V30/V31*V16+V32/AKC2S**2#EZS 
V34= (-SIGN) *2.D0*V14*EP$2/AML)2*V17 
V45=£P$2/AMU2*< EX1*EZ2-EX2*EZ1) 

IF (GAMA2S. LE.O.DO) GC TG 22 

HYPERBOLIC CASE 



V10=GAMA2*D(*N) 

VI 2= 1.00/ (DTANHtVIO )*C-AMA2) 

IF < V10.G6.170.D0 ) GC TC 24 
V 1 3 = D ( M N ) / ( (DSINH(VIO) )**2) 

GO TO 25 

24 V 13 = 0 . 0 0 

25 RE1= (V33+V34)*( V12 + V13 J/2.D0 
ZM1 = -RE 1 

RE3=V45*( V12+V13) 

ZM2=-RE3 
GO TC 23 

TRIGONOMETRIC CASE 

22 V40=D (MN)/( DSIN(GA W A2*C(MN) )**2) 

V4 1=1 . DO / ( DT AN ( G AM A2* C ( MN ) ) *G AM A 2 ) 
RE1=( V33-V34)*( V4C+V41 )/2.D0 

ZM 1= -R = 1 

RE2=-V45*( V40+V41 ) 

ZM2=-RE3 

23 RE11=RE11+ Q E1 
RE2 2=R E 2 2 + RE 3 
ZM 1 1=ZM 1 1+ZM 1 
ZM2 1= ZM2 1 + ZM 2 
GO TC 100 
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19 G1=2.3D6*PI*G 
0MGAM=G1*AMAGS 

gmgao=gi*ho 

OMGAR*GMGAO-CMGAM 
DAMP=Gl*DELTH/( 2.CO*C M EGA ) 

DELT = ( OMGAR**2-OMEGA**2*( l.D0+DAMF**2 ) )**2+ 

X4.D0* (OMGAR*OMEGA*CAMP ) **2 

C CALCULATION OF SUSCEPTIBILITY TEMSQF'S COMPONENTS 

XXXl=OMGAM*OMGAR* ( (OMGAR**2-CMEGA**2 ) + ( CM EGA* DAMP ) ) 
XXX2=OMGAM*QMEGA*CAMP* < CMGAR**2+GM5GA**2* ( 1. DC + CAMP**2 ) ) 
XXZl=OMEGA*OMGAM* (OMGAR**2-OM EGA**2*( 1 . D 0 + 0 AM P** 2 ) ) 

XXZ2 -2. D0*QMEGA** 2*0 MG A N*DAMP*CM GAR 
RE = XXX1*R£1 1 + RE22*XXZ1 
ZM=XXX2*ZMU + ZM21*XXZ2 
904 FORMAT {///) 

WRITE (6,904) 

WRITE (6,904) 

V50=0MEGA*AMU1/(4.CG*T0TPWR*DELT*A) 

ATTEN=-ZM*V50 
D5=8.685889638D-02*ATTEN 
BETAP=BETA+RS*V50 
RAT=BET AP/BET A 
IF (SIGN.EQ.l.DO) GO TO 81 
WRITE (6,898) DB 

393 FORMAT ( 10X , 2 1HBACKW AP C ATTENUAT I CN = , El 1 . 5 , 6H CB/CM//) 
WRITE (6,897) PAT 

897 FORMAT ( 10X , 21H3ET AF/ BETA BACKWARC =,510.4//) 



C S IGN=+1 • DO , PROPAGATION IN FORWARD DIRECTION 
C 



(-Z) 



81 

900 

899 

3001 



S IGN = + 1 • DO 
GG TO 8000 
WRITE (6,900) DB 

FORMAT ( 10X,20HF 3RWARC ATTENUATION^ 
WRITE (6,399) RAT 

FORMAT (1OX,2OHBETAF/0ETA FORWAPQ = 
GO TO 21 



STOP 

END 



E 1 1 • 5 , 6H DB/CM//) 
E 10 • 4/ / ) 
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